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Abstract:  This  paper  identifies  two  issues  with  traditional  root  normal  size  distributions,  which 
are  commonly  fitted  to  experimental  results  for  liquid  atomization  and  sprays.  First,  while  root 
normal  size  distributions  are  typically  expressed  in  terms  of  mass  mean  diameter,  they  do  not 
actually  obtain  the  correct  mass  mean  diameter.  The  error  is  usually  small  but  may  be 
significant  in  some  cases.  Simple  corrective  measures  are  suggested.  Second,  depending  on  the 
chosen  form,  traditional  root  normal  size  distributions  may  exhibit  severe  non  -physical 
singularities.  An  alternative  root  normal  size  distribution  is  suggested  that  experiences  only  mild 
singularities.  A  literature  survey  finds  six  possible  parameter  choices  for  root  normal  size 
distributions.  It  is  found  that  these  parameter  choices  provide  an  adequate  fit  to  a  wide  variety  of 
experimental  data,  except  possibly  for  the  largest  droplets. 

Keywords:  Root  Normal  Size  Distribution,  Simmons  Universal  Size  Distribution,  Spray 
Atomization,  Probability  Density  Function,  Aerosol  Size  Distribution,  Liquid  Fragmentation 

1.  Introduction 

First  introduced  by  Tate  and  Marshall  (1953),  root  nonnal  size  distributions  are  commonly  used 
to  describe  aerosol  size  distributions,  and  appear  in  most  recent  textbooks  on  multiphase  flows, 
atomization,  and  sprays,  e.g.,  Bayvel  and  Orzechowski  (1993),  Liu  (2000),  Brennan  (2005), 
Crowe  (2006),  Ashgriz  (2011). 

Simmons  (1977)  proposed  the  best-known  root  normal  size  distribution.  In  particular,  the 
Simmons  root  nonnal  size  distribution  has  been  featured  in  an  extensive  series  of  papers  by  Dr. 
Gerard  Faeth  and  his  colleagues,  e.g.,  Wu  et.  al.  (1991,  1992  1995);  Ruff  et.  al.  (1992);  Hsiang 
and  Faeth  (1992,  1993);  Wu  and  Faeth  (1993),  Chou  et.  al.  (1997),  Dai  et.  al.  (1998),  Sallam  et. 
al.  (1999,  2006),  Aalberg  et.  al.  (2005),  Lee  et.  al.  (2007),  Miller  et.  al.  (2008). 

As  the  main  advantage  of  root  normal  size  distributions,  the  research  literature  suggests  that 
there  are  only  a  limited  number  of  parameter  choices.  Each  such  parameter  choice  is  sometimes 
referred  to  as  a  universal.  Contrast  this  situation  with  other  distributions,  especially  those  with 
three  or  more  free  parameters,  where  a  small  change  in  the  experimental  data  may  lead  to  a  large 
change  in  the  as-fitted  parameters,  e.g.,  Dumouchel  (2009),  Dumouchel  et.  al.  (2012). 

As  the  main  disadvantage  of  root  normal  size  distributions  -  as  well  as  all  known  alternatives  - 
“no  single  distribution  accurately  fits  even  a  large  fraction  of  the  available  drop  size  data 
research  literature  . . .  [which]  necessitates  trial-and-error  use  of  several  distributions  to  determine 
which  one  best  fits  a  particular  data  set.”  (Ashgriz,  2011).  As  one  main  goal,  this  paper  attempts 
to  clarify  how  close  root  normal  size  distributions,  with  just  a  half-a-dozen  different  parameter 
settings,  can  come  to  fitting  a  wide  variety  of  experimental  data. 
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The  research  literature  offers  few  mathematical  expressions  for  root  nonnal  size  distributions.  As 
a  second  main  goal,  this  paper  considers  two  variants  of  the  root  normal  size  distribution,  one 
common  and  one  rare,  with  complete  mathematical  expressions  given  for  both.  For  example,  this 
treatment  provides  new  algebraic  results  for  the  ratio  of  the  mass  mean  diameter  to  the  Sauter 
mean  diameter,  which  reveal  unexpected  sensitivity  to  an  arbitrary  minimum  droplet  size. 

Root  normal  size  distributions  are  often  expressed  in  normalized  (self-similar)  forms,  i.e.,  with 
the  independent  variable  divided  (and  optionally  the  dependent  variable  multiplied)  by  an 
average  droplet  size.  This  paper  suggests  a  way  to  ensure  that  root  normal  size  distributions 
actually  obtain  the  average  droplet  sizes  implied  by  such  normalized  fonns.  This  requires 
sacrificing  a  free  parameter. 

Finally,  this  paper  confirms  an  earlier  observation,  namely,  that  the  common  variant  of  the  root 
nonnal  size  distribution  experiences  a  severe  singularity  for  small  droplets.  By  contrast,  this 
paper  finds  that  the  rare  variant  has  only  a  mild  singularity.  This  observation  suggests  the  use  of 
the  rare  variant  in  future  work. 

Aerosol  droplets  are  typically  nearly-spherical;  thus  aerosol  size  distributions  are  usually 
expressed  as  functions  of  droplet  diameter.  However,  certain  liquids  tend  to  fonn  highly- 
distended  or  stringy  fragments  including  viscoelastic  liquids  such  as  silicone  oils,  non- 
Newtonian  liquids  such  as  starch  solutions,  and  highly-viscous  liquids  such  as  thick  melts,  e.g., 
Joseph  et.  al.  (1999),  Joseph  et.  al.  (2002),  Theofanous  (2011).  Even  simple  liquids  such  as  water 
may  exhibit  irregularly  shaped  fragments  for  a  brief  period  of  time  immediately  following  an 
atomization  event.  Such  situations  are  best  described  in  terms  of  mass  rather  than  in  terms  of 
diameter.  Thus  this  treatment  gives  expressions  both  for  mass  and  diameter. 

2,  Size  Distributions  and  Transformation  Conditions 

Let  D  be  the  droplet  diameter  and  let  M  be  the  droplet  mass.  Then  there  are  eight  common  ways 
of  expressing  aerosol  size  distributions: 

Fm  (D)  [  Fm  ( M )  ]  is  the  mass  fraction  of  droplets  with  diameters  [masses]  greater  than 
or  equal  to  D  [M]. 

f M  (D)  [  f M  ( M )  ]  is  the  mass  fraction  of  droplets  with  diameters  [masses]  in  a  range  dD 
centered  on  D  divided  by  dD  \dM  centered  on  M  divided  by  dM] 

F(D)  [  F{M )  ]  is  the  number  fraction  of  droplets  with  diameters  [masses]  greater  than  or 
equal  to  D  [M] 

f  ( D )  [  f(  M )  ]  is  the  number  fraction  of  droplets  with  diameters  [masses]  in  a  range  dD 
centered  on  D  divided  by  dD  [dM  centered  on  M  divided  by  dM] 
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This  list  excludes  minor  variations  such  as  1  -  Fu  ,  or  the  use  of  volume  V  instead  of  mass  M. 
The  above  definitions  imply  that  Fu  is  monotone  decreasing  such  that  FM  (0)  =  1  and 
Fm  (go)  =  0 .  In  addition,  fM  is  always  non-negative  such  that: 

oo 

\fM(x)dx  =  1  (la) 

o 

Similarly,  the  above  definitions  imply  that  F  is  monotone  decreasing  such  that  F( 0)  =  1  and 
F(oo)  =  0 .  In  addition,  /  is  always  non-negative  such  that: 

oo 

J/(x)c/x  =  1  (lb) 

0 

In  standard  probability  theory,  F  is  called  a  complementary  cumulative  distribution  function 
(■ CCDF)  and  /  is  called  a  probability  density  function  (PDF). 

As  defined  here,  the  transformation  condition  requires  that  all  eight  forms  given  above  have,  at 
most,  mild  (integrable)  singularities.  It  would  not  make  physical  sense  for  an  aerosol  size 
distribution  to  be  well-behaved  in  one  form  but  to  experience  severe  (non-integrable) 
singularities  in  another  form.  As  seen  below,  traditional  root  normal  size  distributions  do  not 
obey  the  transformation  condition. 

Transforming  between  the  eight  different  forms  given  above  requires  eight  different  equations. 
The  first  of  these  equations  is  as  follows: 


M  =  pCDm 


(2) 


where  p  is  density,  C  is  a  shape  factor,  and  m  is  the  spatial  dimension  (1  <  m  <  3).  For  classic 
aerosols  with  nearly-spherical  droplets,  m  -  3  .  Equation  (2)  may  not  apply  to  certain  highly- 
irregular  liquid  fragments;  this  treatment  specifically  excludes  such  cases. 

The  next  two  transformation  equations  are  as  follows: 


Fm(M)  =  Fm(D ) 

(3) 

F(M)  =  F(D) 

(4) 

which  follow  directly  from  the  definitions  of  the  functions  involved.  The  fourth  transformation 
equation  is  as  follows: 


f  M(x)/(x) 

'  M  W  oo 

J  M(x )  /  ( x)dx 


(5) 
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This  was  first  introduced  by  Brown  (1989)  and  Brown  and  Wohletz  (1995);  see  also  Dumouchel 
(2009).  The  last  four  transfonnation  equations  are  well-known  and  follow  directly  from  the 
definitions  of  the  functions  involved: 


fa AD)  =  -\fM(x)dx', 

D 


FM(M)  =  -\fM(x)dx; 

M 


F(D)  =  -  j  f(x)dx ; 

D 


F(M)  =  -  J  f(x)dx  ; 

M 


(6) 

(7) 

«■>-£ 

(8) 

1 

II 

£ 

(9) 

3.  Average  Droplet  Sizes  and  Self-Similarity  Conditions 

As  noted  above,  average  droplet  sizes  are  often  used  to  nonnalize  size  distribution  functions.  The 
complexity  of  various  mathematical  expressions  depends  on  the  choice  of  average.  Judging  by 
the  complexity  of  the  expressions  given  later  in  Sections  5  to  8,  the  eight  size  distributions 
defined  above  appear  to  be  naturally  associated  with  different  average  sizes.  For  example, 

Fu  ( D )  and  fM  ( D )  are  naturally  associated  with  the  following  two  averages: 


Du„^\Df„,(D)dD 

0 


\Dtf{D)dD 

_0 _ 

00 

\D‘f(D)dD 

0 


\D,f(D)dD 


'^l^FldD  j  i:y  f(D)dn 

n  a 


(10) 


(11) 


which  are  known  as  the  mass  mean  diameter  (MMD)  and  the  Sauter  mean  diameter  ( SMD ), 
respectively.  In  the  research  literature,  it  is  common  to  see  the  ratio  of  these  two  averages: 


R 


M 


F M  avg 

D' 

avg 


(12) 
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Notice  that  D'M  avg  emphasizes  small  droplets  while  DM  avg  emphasizes  large  droplets.  As  a  result, 
1  <  Rm  <  oo  measures  aerosol  size  spread  where  RM  =  1  corresponds  to  the  least  possible  spread 
(i.e.  monodisperse)  while  RM  =  oo  corresponds  to  the  greatest  possible  spread. 

For  another  example,  F(D )  and  f  ( D )  are  naturally  associated  with  the  following  two  averages: 


D.v  =  J  Df(D)dD 

0 


(13) 


D'  = 

avg 


/m(D) 

D 


dD 


(14) 


The  former  is  known  as  the  count  mean  diameter  (CMD).  The  latter  is  non-standard  and 
unnamed.  As  before,  one  can  define  a  ratio  of  these  two  averages: 


R  = 


D. 


avg 


D' 


avg 


(15) 


The  observations  made  earlier  about  Ru  apply  equally  to  R  . 

As  a  third  example,  Fm(M )  and  fM(M )  are  naturally  associated  with  Mm  where: 


\M2f(M)dM 

Mu^=\Mfu(M)dM  =  ^ -  (16) 

0  |  Mf(M)dM 

0 

As  a  fourth  and  final  example,  F(M)  and  f(M)  are  naturally  associated  with  Mavg  where: 


Mu,s=\Mf(M)dM  (17) 

0 


When  transforming  between  the  eight  different  size  distributions  given  in  Section  2,  it  is  useful 
to  have  ratios  such  as  the  following: 


Q 


D 


M  avg 


D 


avg 


(18) 
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(19a) 


SM  = 


M 


M  avg 


pCD 


M  avg 


See  Mavg 


pCD 


avg 


(19b) 


Notice  that  SM  measures  skewness,  where  SM  =  1  if  the  aerosol  size  distribution  is  evenly 
balanced  between  small  and  large  droplets,  e.g.,  a  uniform  size  distribution.  Similar  observations 
apply  to  S. 

In  the  simplest  possible  case,  self-preserving  or  self-similar  aerosol  size  distributions  can  be 
written  in  terms  of  the  following  ratio: 


M 

?7  = - 

M 

ly±  M  avg 

This  was  first  observed  by  Friedlander  and  Wang  (1966)  for  populations  of  droplets  that  collide 
with  each  other;  see  also  Lee  et.  al.  (1984),  Spicer  and  Pratsinis  (1996),  and  Lehtinen  and 
Zachariah  (2001).  For  example,  Spicer  and  Pratsinis  (1996)  give  the  following  definition:  “When 
the  steady-state  size  distributions  scaled  by  the  average  particle  volume  [or  mass]  collapse  onto  a 
single  size  distribution,  this  distribution  is  tenned  self-preserving.” 

Using  the  transfonnations  given  in  Section  2,  self-similar  aerosol  size  distributions  can  equally 
well  be  written  in  terms  of  the  following  ratio: 


In  fact,  as  already  noted,  root  nonnal  size  distributions  are  commonly  written  in  such  simple  self¬ 
preserving  (nonnalized)  forms.  An  obvious  condition  is  that,  if  an  aerosol  size  distribution  is 
written  in  terms  of  a  given  average,  it  should  actually  obtain  that  average.  However,  as  seen 
below,  traditional  root  normal  size  distributions  do  not  obey  this  condition. 

4.  Modifications  for  Minimum  Droplet  Sizes 

The  vast  majority  of  experiments  impose  a  minimum  droplet  size  Z)min  or,  equivalently, 

Mniin  =  pCD'f  due  to  optical  resolution  or  other  practical  factors.  If  the  minimum  droplet  size 
is  large  enough,  the  aerosol  size  distribution  should  be  modified  accordingly.  For  example,  FM 
should  be  replaced  by  FM  where  FM  (Dmin )  =  1  and  FM  (go)  =  0 .  This  is  true  if: 


Fm(D)  = 


Fm{D ) 

FxfDmm) 


(20) 
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Similarly,  fM  should  be  replaced  by  fM  where: 


]fM(D)dD  =  1 


This  is  true  if 


/m(A>  = 


fM(P) 

Fm  (A™  ) 


(21) 


As  another  example,  DMavg  should  be  replaced  by  PMavs  where: 


M  avg 


OO  |  CO 

=  I  D~fu  (D)dD  =  — — —  |  Df„  (D)dD 


(22) 


Similarly,  D'l/ais  should  be  replaced  by  P'Mavg  where: 


M  avg 


D[, 


FM  (Pmm  ) 


M  avg  oo 


J  j 


D 


fu(D) 

D 


dD 


(23) 


As  a  final  example,  RM  should  be  replaced  by  RM  where: 


D  _  D M  avg 

"  = 

avg 


(24a) 


or: 


Rm  ~ 


1 


FM  (Pnan  )  n 


CO  00 

\DfM(D)dD  J  -fu(D)dD 


(24b) 


Alternatively,  when  DM  is  finite  while  D'M  is  not,  it  is  common  to  see  mixed  expressions 
such  as  the  following: 


-  A 
Rm  -  = 


D* 

^ M  avg 


(25) 
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5.  Traditional  Root  Normal  Size  Distributions  (Type  I) 

The  traditional  root  normal  size  distribution  may  be  written  as  follows: 


Table  1  recasts  this  size  distribution  into  eight  equivalent  forms  using  the  expressions  given  in 
Section  2.  Unfortunately,  there  are  no  known  analytical  expressions  for  two  of  these  fonns  - 
F(D )  and  F(M)  -  meaning  they  must  be  detennined  numerically. 


Table  1  a.  Root  normal  distributions  of  Type  I  expressed  in  terms  of  FM  and  fM 


Table  lb.  Root  normal  distributions  of  Type  I  expressed  in  terms  of  F  and /. 


As  described  in  the  next  section,  the  traditional  root  normal  size  distribution  suffers  from  a  mild 
singularity  in  fM  but  a  severe  singularity  in  /  at  the  origin.  This  singularity  has  been  previously 
observed  by,  for  example,  Babinsky  and  Sojka  (2002),  who  called  it  “a  gradient  catastrophe  near 
zero.” 
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It  is  convenient  to  express  the  properties  of  the  root  nonnal  size  distribution  in  terms  of  the 
following  integral: 


1  r  M/2  (Jx-a) 

-  x’  exp  - - ^ 

, _ 2  J  ^  Irr2 


Ln<j  0 


There  are  analytic  expressions  for  Wi  if  and  only  if  i  is  a  non-negative  integer;  see  Table  2. 


Table  2.  Six  examples  of  integrals  defined  by  Equation  (27)  where  erfc  is  the  complementary  error  function. 


9 


By  Equation  (1): 


(28) 

(29) 


A  =  W0 
B  =  Wm 


By  Equations  (12),  (15),  and  (18): 


R 


M 


w  xwx 

~w 


(30) 


R  = 


w  w 

vv  -m+Vr  -m- 1 


(31) 


Q  = 


K  wm 
K  wm+l 


(32) 


Because  of  the  singularity  at  the  origin,  Equations  (29)  to  (32)  require  imposing  a  minimum 
droplet  size  Z)niin  .  The  minimum  droplet  size  must  be  applied  to  those  Wt  with  a  negative  index 

i  and,  optionally,  to  all  Wt .  In  other  words,  B  should  be  replaced  by  B  ,  RM  should  be  replaced 
by  Rm  ,  R  should  be  replaced  by  R  ,  and  Q  should  be  replaced  by  Q  as  described  in  Section  4. 


Figure  1  shows  how  RM  depends  on  Dmin  and  a.  Notice  that  RM  is  relatively  insensitive  to  the 
arbitrary  choice  of  Dmm  for  0  <  a  <  0.2 .  However,  RM  may  be  extremely  sensitive  to  Z)nim  for 
larger  a .  It  is  a  common  practice  in  the  research  literature  to  specify  Ru  rather  than  <7  . 
However,  this  may  be  somewhat  confusing  given  that  RM  depends  on  an  arbitrary  Dmm  ,  at  least 
for  a  >  0.2  For  most  cases  seen  in  the  research  literature,  Dmin  /DMavg  is  somewhere  between 
0.1  and  0.2.  Figure  2  shows  that  RM  may  vary  substantially  over  this  range  when  a  >  0.2  . 
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1.5 


Omin/^Movg-O-Ol  / 


0.1/ 


Figure  1.  The  ratio  RM  for  traditional  (Type  I)  root  normal  size  distributions  with  a= 1. 

If  Fm(D )  and  fM(D)  are  self-similar,  then  DM  in  Table  la  must  agree  with  Du  as 
defined  by  Equation  (10).  This  is  true  if: 

(33) 


According  to  Table  2,  this  can  be  written  as  follows: 


2  2  2 era  exp(-a2 /2cr2)  , 

<j-  +a  +^=—^ - =—^  =  1 

V2 n  erfc(-a/ v2<r) 


Equations  (33)  and  (34)  have  not  appeared  in  the  research  literature  before.  Rather,  past 
treatments  assumed  a  - 1 .  Figure  2  compares  Equations  (33),  Equations  (34),  and  the  traditional 
choice  of  a  - 1 . 
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0.96 


o  0.92 
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0  0.1  0.2  0.3  0.4  0.5 


a 


Figure  2.  The  relationship  between  parameters  a  and  a  required  to  ensure  that  traditional  (Type  I)  root 
normal  size  distributions  have  the  correct  count  mean  diameter.  In  the  legend,  ‘self  similar’  refers  to 
Equation  (33)  while  ‘approximate’  refers  to  Equation  (34). 

Figure  3  illustrates  the  effects  of  changing  a.  More  specifically,  Figure  3  compares  the  effects  of 
using  a  fixed  value  of  a  =  1  versus  Equation  (33).  In  this  example,  the  choice  a  - 1  results  in 
mass  mean  diameter  that  is  about  6%  too  large  while  the  choice  a  -  0.9708  results  in  the  correct 
mass  mean  diameter. 


D/&M  avg 

Figure  3.  Traditional  (Type  I)  root  normal  size  distributions  with  o=  0.24  and  a= 1  vs.  «=0.9708.  The  latter 
value  of  a  was  chosen  using  Equation  (33)  to  ensure  the  correct  mass  mean  diameter. 
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If  Fm(M )  and  fM(M )  are  self-similar,  then  Mm  in  Table  la  must  agree  with  Mm  as 
defined  by  Equation  (19).  This  is  true  if: 


W 

o  _  m 

M  =  K 

For  example,  if  m  =  3  and  Equation  (33)  is  exactly  true,  then  Equation  (35)  becomes: 

SM  =  a4  +\4a2 a2  +  33a4  -2a2 a4  -18a6 

6.  Variant  of  Traditional  Root  Normal  Size  Distributions  (Type  IB) 

Suppose  the  traditional  (Type  I)  root  nonnal  size  distribution  is  modified  as  follows: 


fu(D)  = 


ADavg  V 2ko 2  ]D!Da 


r  exp 


-a)2 


(35) 


(36) 


Notice  that  Equation  (36)  is  the  same  as  Equation  (26),  except  that  the  mass  mean  diameter 
DMavg  has  been  replaced  by  the  count  mean  diameter  D  .  Equation  (36)  can  be  recast  in  eight 

different  fonns.  The  results  are  identical  to  those  given  earlier  in  Table  1  after  replacing  DMavg 
by  Davg,  MMavg  by  Mavg ,  and  SM  by  S. 

Most  properties  of  Equation  (36)  are  exactly  the  same  as  those  of  Equation  (26).  In  particular, 
Equations  (28)  to  (32)  are  the  same  as  before.  Only  the  self-similarity  conditions  change  to 
become: 


W_m  =  W_m+1  (37) 

W 

S  =  -5-  (38) 

W 

-m 

Unlike  Equations  (33)  and  (35),  Equations  (37)  and  (38)  require  imposing  a  minimum  droplet 
size  Dmm  .  This  extra  complication  makes  Type  IB  root  normal  size  distributions  unattractive. 
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7.  Alternative  Root  Normal  Size  Distributions  (Type  II) 


Consider  the  following  alternative  root  normal  size  distribution: 

1  1  1 


f(D)  = 


AD Mavg  yh/ra2  -^D/D 


rCXp 


M  avg 


DiQDID^ -ay 


(39) 


Notice  that  Equation  (39)  is  identical  to  Equation  (26)  except  that  fM  has  been  replaced  by  /  . 
While  this  alternative  root  normal  distribution  has  been  mentioned  in  the  research  literature 
before,  e.g.,  Crowe  (2006),  its  properties  have  never  been  well-described.  Table  3  recasts  this 
size  distribution  into  eight  equivalent  forms  using  the  expressions  given  in  Section  2. 

As  specific  examples,  Table  4  gives  expressions  for  FM  and  fM  for  m  =  1 ,  2,  and  3.  There  are 
analytical  expressions  for  FM  if  and  only  if  m  is  an  integer. 

As  seen  in  Figure  4a,  when  expressed  in  terms  of  fM  ,  the  traditional  root  nonnal  size 
distribution  has  a  mild  singularity  at  the  origin  while  the  alternative  root  nonnal  size  distribution 
has  no  singularity.  As  seen  in  Figure  4b,  when  expressed  in  tenns  of  /  ,  the  traditional  root 
nonnal  size  distribution  has  a  severe  singularity  at  the  origin  while  the  alternative  root  normal 
size  distribution  has  only  a  mild  singularity.  In  other  words,  the  alternative  root  normal  size 
distribution  satisfies  the  transfonnation  condition  proposed  in  Section  2  while  the  traditional  root 
nonnal  size  distribution  does  not. 
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Table  3a.  Root  normal  distributions  of  Type  IIA  expressed  in  terms  of  F  and /. 


Table  4a.  Root  normal  distributions  of  Type  IIA  expressed  in  terms  of  FM  and  jM  for  m  =  1 

T7  rTW  1/2  2\  f  ^ ^Mavg  a 

FM(D)  =  ~(a  +a  )erfC  - 7= - 


Z7  { A  J\  1  (  2  ,  2  ^  f  mM  I M  Mavg  a 

Fu (M)  =  -(a  +o  )erfc  - r - 

B  v  2o 


+vS^^7^+a 


—QD/DMaig  ~ay 


+  ^b^SmM/Mm°V8  +  a)exp 


1  1  rpyp; -  QD/DMavg-aY  sm  1  -  W^/M 

DIDMavs^ - pp -  fM(M)  =  ~ - - ^^JSMM/MMavexp - — 


FFfMavg  \27T<72 


Mavg  V2 no1 


Table  4b.  Root  normal  distributions  of  Type  IIA  expressed  in  terms  of  FM  and  fM  for  m  =  2 


1  ,  4  ,  £  2  2  ,  (D'DMavg)  ~a 


FM(D)  = —(a  +6a  a  +3cr  )erfc 


2a  \  D  2  £>  2  0  2J  D 

h — .  —  -  +  Q - h  (q  +  3cr  )  - 

\27rB  F  M  av2  >  F M  av2  .  F M  a\ 


1  /  4  .  /  2_2  .  p  ^ M  avg 


Fm(M)  =  —(a  +  6a  cr  +3cr  )erfc 


(■ SMM/MMavgy 

Jla 


+ ~^L~  M-  4+JM_  2+(fl2+3(72 

yinB  M  Mavs  M  Mav2  M  A 


+  a 3  +5 acr1  I  expl  - 


((D  /  DM  y12  -  a)2 


M  avg  7 


+  a3  +5acr2  I  expl  - 


{(SMM/MMYl4-a)2 


Table  4c.  Root  normal  distributions  of  Type  IIA  expressed  in  terms  of  FM 


1 


FM{D)  =  -{a  +\5aaz  +  45<Tct4  +  15a6) 
B 


x  erfcs 


1 


44c 


f  Y/2 

D 


\^Mavg  J 


-a 


2(7 


444b 


D 


\DMavg  j 


+  u 


D 


\^Mavg  J 


+  (a2  +  5cr2) 


f  D 


D 


+  (a3  +  9aa2)-  '  '  ~4  '  1  ° ~2 -2  .  1 


D 


■  +  (a4  +12a2<j2  +  15cr4) 


\J-^Mavg  J 


r  D 


M  avg 


y^Mavg  J 


+  a3  +14a3cr1  +33  aa 


exp 


2cr2 


D 


y^Mavg  J 


-a 


f M  (P)  = 


BDavg  44 


KG 


D 


y^Mavg  J 


exp 


2a2 


D 


y  Dm  avg  ) 


-a 


for  m  =  3 


Fm(M)  =  —(a6  +  15a4cr2  +  45a2a4  +  15<x6) 


xerfc  >  \^L 

V  2cr  M. 


\2nB  M  M  M  Mav2  M 


c  /If  l3  c 

+  (a 3  +9aa2)  — —  +  (a4  +12 a2a2  +15 a4)  - 

l  M Mavg  J  l  ^ 


1  C  M  I6 

+  a5  +I4a3a2  +33acr4  expJ - ^  — — 


0.0000001  0.000001  0.00001  0.0001  0.001  0.01  0.1  1  10 


D/Dm  avg 

(a.) / m(D )  for  Type  I  vs.  Type  II  root  normal  size  distributions  in  a  linear-log  plane. 


0.01  0.1  1  10 

D/Dm  avg 

(b.)  f(D)  for  Type  I  vs.  Type  II  root  normal  size  distributions  in  a  linear-log  plane. 

Figure  4.  The  Type  I  root  normal  size  distribution  experiences  a  mild  singularity  in /m(D)  and  a  severe 
singularity  in  f(D).  The  Type  II  root  normal  size  distribution  does  not  experience  any  singularity  in  f  m(D ) 
and  experiences  only  a  mild  singularity  in  f(D).  The  parameter  settings  for  the  Type  I  distribution  are 
o=0.238  and  a= 1.  The  parameter  settings  for  the  Type  II  distribution  are  o=0.27,  «=0.5019,  and  m=  3  where 
the  value  of  a  is  chosen  to  ensure  the  correct  mass  mean  diameter. 
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By  Equation  (1): 


By  Equations  (12),  (15),  and  (18): 


A  =  W0 
B  =  Wm 


w  w 

VV  m-V'  n 


m+ 1 


W, 


R  =  W 


w, 


Q  = 


Wo  wm+l 


(40) 

(41) 

(42) 

(43) 

(44) 


Equation  (43)  requires  imposing  a  minimum  droplet  size  Dmm  .  This  minimum  droplet  size 
should  be  applied  to  W_x  and,  optionally,  W0  and  If] .  In  other  words,  R  should  be  replaced  by  R 
as  described  in  Section  4. 


Figure  5  shows  RM  as  a  function  of  a.  Unlike  Type  I  distributions,  Type  II  distributions  have  a 
one-to-one  analytical  relationship  between  Ru  and  a  .  Notice  that  the  Type  II  results  agree  well 
with  the  Type  I  results  when  the  latter  has  Dmin  / DMavg  =  0. 1 5  . 
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Figure  5.  The  ratio  RM  for  alternative  (Type  II)  root  normal  size  distributions  with  a  chosen  to  ensure  self¬ 
similarity.  For  comparison,  results  for  traditional  (Type  I)  root  normal  distributions  with  a=l  are  also  shown. 

If  Fm(D )  and  fM(D )  are  self-similar,  then  DM  in  Table  3  must  agree  with  DM  as  defined 
by  Equation  (10).  This  is  true  if: 

Wm  =  Wm+1  (45) 

Figure  6  shows  how  a  varies  with  o  per  Equation  (45). 
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1 


a 


Figure  6.  The  relationship  between  parameters  a  and  a  required  to  ensure  that  Type  II  root  normal  size 
distributions  have  the  correct  mass  mean  diameter.  This  assumes  that  m= 3.  This  is  a  plot  of  Equation  (45). 

If  Fm(M )  and  fM{M )  are  self-similar,  then  MMavg  in  Table  3  must  agree  with  MMavg  as 
defined  by  Equation  (19).  This  is  true  if: 


W 

c  — 

^  M  ~ 


W„, 


(46) 


8.  Variant  of  Alternative  Root  Normal  Size  Distributions  (Type  IIB) 

Suppose  the  alternative  (Type  II)  root  nonnal  size  distribution  is  modified  as  follows: 


f(D)  = 


1 


1 


1 


ADavg  a/ 2ttct2  p/ Davg 


rcxp 


2  cr 


-aY 


(47) 


Notice  that  Equation  (47)  is  the  same  as  Equation  (39),  except  that  the  mass  mean  diameter 
DMavg  has  been  replaced  by  the  count  mean  diameter  D  .  Equation  (47)  can  be  recast  in  the 

usual  eight  forms.  The  results  are  identical  to  those  given  earlier  in  Table  3  after  replacing  DMavg 
by  Davg,  MMavg  by  Mavg,  and  SM  by  S. 

Most  properties  of  Equation  (47)  are  exactly  the  same  as  those  of  Equation  (39).  In  particular, 
Equations  (40)  to  (44)  are  the  same  as  before.  Only  the  self-similarity  conditions  change  to 
become: 
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(48) 


^0=^1 

W 

S  —  — — 
^0 


(49) 


Notice  that  Equations  (48)  and  (49)  for  Type  IIB  distributions  are  the  same  as  Equations  (33)  and 
(35)  for  Type  I  distributions,  after  replacing  SM  by  S.  Thus  the  earlier  discussion  applies  here, 
including  Figures  2  and  3.  . 

Figure  7  gives  an  example  showing  that,  with  the  right  parameter  choice,  Type  II  and  IIB  root 
nonnal  size  distributions  are  nearly  identical.  Notice  that  the  Type  IIB  distribution  is  plotted  as 
QfM  ( QD )  .  As  their  main  advantage,  Type  IIB  distributions  have  somewhat  simpler  self- 
similarity  conditions  than  Type  II  distributions. 


0  0.5  1  1.5  2  2.5  3 

D/Dm  avg 

Figure  7.  An  example  of  Type  II  vs.  Type  IIB  root  normal  size  distributions  for  m=3.  The  parameter  a  in  both 
cases  is  chosen  to  ensure  self-similarity,  i.e.,  the  Type  II  has  the  correct  MMD  and  the  Type  IIB  has  the 
correct  CMD. 

9.  Experimental  Evidence 

Table  5  lists  six  root  nonnal  size  distributions  found  in  a  literature  search.  Each  reference  given 
in  Table  5  typically  provides  three  root  normal  size  distributions  —  an  upper  bound,  a  lower 
bound,  and  an  average  -  to  describe  anywhere  between  four  and  2,000  separate  tests.  For  the 
references  marked  by  an  asterisks,  the  given  root  nonnal  size  distribution  is  a  lower  bound  to 
multiple  sets  of  test  data.  For  the  remainder,  the  given  root  normal  size  distribution  is  an  average 
of  multiple  sets  of  test  data. 
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Table  5.  Six  examples  of  traditional  (Type  I)  root  normal  size  distributions  found  in  the  research  literature. 
The  values  given  in  bold  are  specified  in  the  original  reference.  The  other  values  are  derived.  All  root  normals 
are  averages  of  multiple  sets  of  test  data  except  where  noted. _ 


a 

a 

Rm 

Reference(s) 

1 

0.081 

1 

1.02 

Chou  &Faeth  (1998)  (*) 

2 

0.11 

1 

1.04 

Sallam  et.  al.  (2006)  (*) 

3 

0.17 

1 

1.10 

Wu  et.  al.  (1991)  (*);  Sallam  et.  al.  (2006) 

4 

0.20 

1 

1.15 

Empie  et.  al.  (1993,  1995,  1997) 

5 

0.238 

1 

1.20 

Simmons  (1977);  Wu  et.  al.  (1991);  Ruff  et.  al.  (1992);  etc. 

6 

0.286 

1 

1.5  (t) 

Speilbauer  et.  al.  (1989);  Adams  et.  al.  (1990);  Loebker  &  Empie  (2001) 

(f)  Highly  sensitive  to  the  choice  of  minimum  droplet  size. 
(*)  Lower  bound  to  multiple  sets  of  test  data  in  terms  of  RM  . 


In  general,  the  curve  fits  in  Table  5  were  derived  using  a  traditional  procedure  that  involves,  first, 
transforming  the  test  data  into  a  plane  in  which  the  presumed  size  distribution  is  linear  and, 
second,  using  least  squares  to  estimate  the  slope.  In  the  case  of  power  laws,  Clauset  et.  al. 

(2009)  has  criticized  this  approach  as  follows:  “ commonly  used  methods  for  analyzing  power-law 
data,  such  as  least-squares  fitting,  can  produce  substantially  inaccurate  estimates  of  parameters 
...  Even  in  cases  where  such  methods  return  accurate  answers  they  are  still  unsatisfactory 
because  they  give  no  indication  of  whether  the  data  obey  a  power  law  at  a//.”  The  same  criticism 
applies  in  general,  including  to  root  nonnal  size  distributions.  However,  regardless  of  such 
concerns,  some  of  the  root  normal  fits  listed  in  Table  5  -  especially  Simmons  universal  -  have 
been  successfully  applied  to  a  wide  variety  of  tests  for  decades. 

Pimentel  et.  al.  (2010)  compiled  42  test  results  for  spray  atomization,  including  both  their  own 
original  results  and  those  found  in  a  literature  survey  done  by  Paloposki  (1994).  This  is  a  neutral 
source  of  test  data,  in  the  sense  that  most  of  it  has  not  been  previously  fitted  by  root  nonnal  size 
distributions.  As  shown  in  Table  6,  all  of  these  test  results  can  be  reasonably  well-fitted  by  one 
of  the  six  root  nonnal  size  distributions  listed  in  Table  5,  with  the  exception  of  three  tests  due  to 
Bhatia  et.  al.  (1988),  which  are  omitted. 

The  size  distributions  given  in  Table  6  were  chosen  to  minimize  the  average  enor  in  four 
different  views,  namely,  fM  (D)  in  the  linear-linear  plane,  fM  (D)  in  the  log-log  plane,  f  (D)  in 
the  linear-linear  plane,  and  f  ( D )  in  the  log-log  plane.  In  essence,  each  view  implies  a  different 
weighting  on  the  error.  In  most  cases,  a  maximum  droplet  size  was  imposed  prior  to  choosing  the 
size  distribution.  This  affects  quantities  such  as  the  mass  mean  diameter.  The  maximum  droplet 
size  addresses  three  issues:  first,  it  balances  out  the  effects  of  minimum  droplet  sizes  inherent  to 
most  tests;  second,  it  avoids  random  variations  caused  by  small  samples  of  large  droplets  seen  in 
some  tests;  and,  third,  it  avoids  divergence  from  root  normal  size  distributions  for  the  largest 
droplets  seen  in  some  tests. 
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Table  6.  Experimental  results  for  spray  atomization  given  by  Pimentel  et.  al.  (2010).  As  indicated  in  the  first 
and  second  columns,  all  results  can  be  reasonably  well-fitted  by  one  of  the  root  normal  size  distributions 
listed  in  Table  5,  except  for  those  from  Bhatia  et.  al .  (1988),  which  are  omitted. _ _ _ 


a 

Reference(s) 

Liquid 

Spray  System 

No. 

% 

1 

0.081 

Pimentel  et.  al.  (2010) 

Jet  Fuel 

‘Delavan’ 

3 

7.8 

2 

0.11 

Pimentel  et.  al.  (2010) 

Jet  Fuel 

‘Bosch’ 

4 

35.9 

ft 

ft 

‘BETE’ 

9 

ft 

ft 

‘Delavan’ 

1 

3 

0.17 

Pimentel  et.  al.  (2010) 

(unknown) 

‘LaVision’ 

3 

35.9 

Paloposki  &  Fagerhohn  (1986) 

Fuel  Oil 

Hollow  Cone 

5 

Li  &  Tankin  (1987) 

Water 

Solid  Cone 

4 

Tishkoff  (1979) 

Water 

Hollow  Cone 

2 

4 

0.20 

Turner  &  Moulton  (1953) 

P-naphthol 

Swirl  Jet 

4 

10.2 

5 

0.238 

Tate  &  Olson  (1962) 

Water 

Solid,  Hollow  Cone 

3 

10.2 

Houghton  (1941) 

(unknown) 

(unknown) 

1 

6 

0.286 

- 

- 

- 

0 

0 

The  following  subsections  review  the  experimental  evidence  found  in  the  research  literature  for 
each  of  the  size  distributions  listed  in  Table  5.  In  addition,  the  following  subsections  compare 
test  data  to  the  chosen  root  normal  size  distributions  for  each  case  listed  in  Table  6. 

9.1  Experimental  Evidence  for  o  =  0.081 

As  an  example,  Chou  and  Faeth  (1998)  studied  aerodynamic  breakup  of  water  and  glycerol- 
water  droplets.  They  found  that  bag-type  breakup,  excluding  the  basal  ring,  led  to  Type  I  root 
nonnal  size  distributions  with  a  lower  bound  of  RM  ~  1.02 ,  an  upper  bound  of  RM  «  1.10,  and 
an  average  of  RM  ~  1 .04  .  The  two  tests  with  the  most  glycerol  (42%,  63%)  did  not  consistently 
obtain  the  lower  bound,  upper  bound,  or  average.  However,  the  two  tests  with  the  least  glycerol 
(0%,  21%)  approximately  obtained  the  lower  bound  of  RM  »  1.02  or,  equivalently,  a  ~  0.081 
and  a  = 1 . 

As  another  example,  Figure  8  compares  Type  I  and  II  root  normal  distributions  against 
experimental  data  obtained  by  Pimentel  et.  al.  (2010)  using  several  different  Delavan  atomizers. 
The  Type  I  and  II  distributions  are  essentially  identical.  The  root  normal  distributions  only  agree 
with  the  experimental  results  for  the  smallest  droplets.  The  larger  droplets  appear  to 
approximately  obey  a  power  law  rather  than  a  root  nonnal  distribution,  as  indicated  by  linearity 
in  the  log-log  plane. 
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(a.)  Log-log  plane 


D/Dm  avg 


D/D, M  avg 


(b.)  Linear-linear  plane 


Figure  8.  Type  I  vs.  Type  II  root  normal  distributions  compared  to  experimental  results  obtained  by  Pimentel 
et.  al.  (2010)  using  Delavan  atomizers. 
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9.2  Experimental  Evidence  for  o  =  0.11 

As  an  example,  Sallam  et.  al.  (2006)  studied  primary  breakup  of  round  aerated  jets  of  water, 
ethyl  alcohol,  and  79%  glycerol  solutions  in  supersonic  crossflows.  They  found  Type  I  root 
nonnal  size  distributions  with  a  lower  bound  of  RM  ~  1.04  and  an  upper  bound  of  RM  -1.10. 
Two  of  the  five  tests  were  somewhere  between  the  upper  and  lower  bounds.  One  of  the  five  tests 
approximately  obtained  the  upper  bound.  Two  of  the  five  tests  approximately  obtained  the  lower 
bound;  more  specifically,  water  and  glycerol  with  jet  diameters  of  0.5mm  approximately 
obtained  RM  «  1.04  or,  equivalently,  a  ~  0.1 1  and  a  -  1 . 

As  another  example,  Figure  9  compares  Type  I  and  II  root  normal  distributions  against 
experimental  data  obtained  by  Pimentel  et.  al.  (2010)  using  a  common  rail  Bosch  system  and  a 
BETE  XA-PR200  nozzle.  The  Type  I  and  II  distributions  are  essentially  identical.  The  Type  I 
and  II  distributions  agree  reasonably  well  with  the  experimental  results,  except  for  the  largest 
droplets.  As  before,  in  the  experiments,  the  largest  droplets  appear  to  approximately  obey  a 
power  law  rather  than  a  root  normal  distribution. 

9.3  Experimental  Evidence  for  a  =  0.17 

As  an  example,  Wu  et.  al.  (1991)  studied  primary  breakup  of  round  turbulent  and  non-turbulent 
jets  of  water,  glycerol- water  solutions,  and  n-Heptane  in  still  air.  They  found  Type  I  root  nonnal 
size  distributions  with  a  lower  bound  of  RM  -1.1,  an  upper  bound  of  RM  «  1.5 ,  and  an  average 
of  Rm  -1.2.  Three  of  the  eighteen  tests,  those  with  relatively  low  jet  speeds,  approximately 
obtained  the  lower  bound  of  RM  « 1.1  or,  equivalently,  a  ~  0.17  and  a  =  1 . 

As  another  example,  Figure  10  compares  Type  I  and  II  root  normal  distributions  against 
experimental  data  obtained  by  Li  and  Tankin  (1987)  and  Tischkoff  (1979).  The  Type  I  and  II 
distributions  agree  well  except  for  the  smallest  droplets.  Where  the  two  distributions  disagree, 
the  experimental  results  slightly  favor  the  Type  II  distribution. 

9.4  Experimental  Evidence  for  a  =  0.20 

As  an  example,  Empie  et.  al.  (1993,  1995,  1997)  conducted  over  140  different  experiments 
involving  black  liquor  sprayed  through  vibrating  and  stationary  splashplate,  V-jet,  and  swirl  cone 
nozzles.  They  found  Type  I  root  nonnal  size  distributions  with  a  lower  bound  of  cr  =  0. 15 ,  an 
upper  bound  of  cr  =  0.28 ,  and  an  arithmetic  average  of  cr  =  0.20  .  The  tables  and  figures  in 
Empie  et.  al.  (1993)  indicate  that  at  least  16  of  the  experiments  obtained  cr  =  0.20  almost 
exactly. 

As  another  example,  Figure  1 1  compares  Type  I  and  II  root  normal  distributions  against 
experimental  data  obtained  by  Turner  and  Moulton  (1953).  The  Type  I  and  II  distributions  agree 
well  except  for  the  smallest  droplets.  Where  the  two  distributions  disagree,  the  experimental 
results  mostly  favor  the  Type  II  distribution. 
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Figure  9.  Type  I  vs.  Type  II  root  normal  distributions  compared  to  experimental  results  obtained  by  Pimentel 
et.  al.  (2010)  using  a  common  rail  Bosch  system  and  a  BETE  XA-PR200  nozzle. 
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Figure  10.  Type  I  vs.  Type  II  root  normal  distributions  compared  to  experimental  results  obtained  by  Li  and 
Tankin  (1987)  and  Tischkoff  (1979).  Both  sets  of  experimental  results  are  as  reported  by  Pimentel  et.  al. 
(2010). 
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(a.)  Log-log  plane 


D/DM  avg 


(b.)  Linear-linear  plane 


D/Dm  avg 


Figure  11.  Type  I  vs.  Type  II  root  normal  distributions  compared  to  experimental  results  obtained  by  Turner 
and  Moulton  (1953).  The  experimental  results  are  as  reported  by  Pimentel  et.  al.  (2010). 
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9.5  Experimental  Evidence  for  a  =  0.238 

As  a  well-known  example,  Simmons  (1977)  examined  “over  2,000  separate  tests  on  about  100 
different  nozzle  designs.”  According  to  Spielbauer  et.  al.  (1989),  the  entire  set  of  tests  yielded 
Type  I  root  nonnal  sized  distributions  with  an  average  cr  of  about  “0.24,  with  a  maximum  17% 
greater,  and  a  minimum  20%  smaller.”  Assuming  a  -  1 ,  this  roughly  corresponds  to  a  lower 
bound  of  Rm  «  1.1 ,  an  upper  bound  of  RM  «  1.5 ,  and  an  average  of  RM  «  1.2  .  For  a  subset  of 
200  tests  “selected  almost  at  random,”  Simmons  (1977)  reported  a  much  narrower  range  of 
Rm  ~  1.20  ±  5% .  Simmons  (1977)  notes  that  RM  =1.20  corresponds  to  cr  =  0.238  and  a  - 1 . 

As  another  example,  Ruff  et.  al.  (1992)  studied  primary  and  secondary  breakup  of  laminar  and 
turbulent  round  water  jets  in  still  air.  They  found  Type  I  root  normal  size  distributions  with  a 
lower  bound  of  RM  «  1 . 1 ,  an  upper  bound  of  RM  «  1 .5  ,  and  an  average  of  RM  «  1 .2  .  For  both 
laminar  and  turbulent  flows,  they  found  RM  ~  1 .2  at  the  jet  surface.  For  laminar  flows,  they 
found  Ru  increased  moving  away  from  the  jet  surface,  e.g.,  1.10,  1.14,  1.22,  1.24.  For  turbulent 
flows,  they  found  RM  decreased  moving  away  from  the  jet  surface,  e.g.,  1.27,  1.22,  1.18,  1.19. 

As  a  third  example,  Hsiang  and  Faeth  (1992,  1993)  and  Chou  et.  al.  (1997)  studied  stripping-type 
breakup  of  liquid  droplets  moving  through  still  air.  In  their  schematic,  aerodynamic  forces  cause 
a  boundary  layer  to  form  on  the  surface  of  a  droplet.  Once  the  boundary  layer  becomes  thick 
enough,  it  gradually  sheds  children  droplets,  whose  mass  mean  diameters  grow  with  the  square 
root  of  time.  Hsiang  and  Faeth  (1992,  1993)  represented  all  children  from  a  given  parent  by  a 
single  size  distribution,  while  Chou  et.  al.  (1997)  represented  the  children  from  a  given  parent  by 
either  three  or  four  different  size  distributions.  All  told,  Hsiang  and  Faeth  (1992,  1993)  and  Chou 
et.  al.  (1997)  reported  stripping  breakup  results  for  fifteen  different  for  combinations  of  liquid 
types  and  initial  relative  droplet  velocities.  Overall,  they  found  Type  I  root  normal  size 
distributions  with  a  lower  bound  of  RM  —1.1,  an  upper  bound  of  RM  «  1.5 ,  and  an  average  of 
Rm  ~  1 .2  .  The  test  results  were  densely  scattered  between  the  upper  and  lower  bounds  without 
any  obvious  clustering  around  the  boundaries  or  the  average. 

As  a  final  example,  Figure  12  compares  Type  I  and  II  root  nonnal  distributions  against 
experimental  data  obtained  by  Houghton  (1941)  and  Tate  and  Olson  (1979).  The  Type  I  and  II 
distributions  agree  well  except  for  the  smallest  droplets.  Where  the  two  distributions  disagree, 
the  experimental  results  clearly  favor  the  Type  II  distribution.  In  the  one  result  from  Houghton 
(1941),  the  largest  droplets  appear  to  approximately  obey  a  power  law  rather  than  a  root  normal 
size  distribution. 
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Figure  12.  Type  I  vs.  Type  II  root  normal  distributions  compared  to  experimental  results  obtained  by 
Houghton  (1941)  and  Tate  and  Olson  (1962).  The  experimental  results  are  as  reported  by  Pimentel  et.  al. 
(2010). 


31 


As  noted  earlier,  past  researchers  generally  fitted  experimental  results  in  the  probability-square- 
root  plane,  i.e.,  in  a  plane  in  which  1  -  FM  ( D )  is  linear.  To  check  the  consistency  between  our 
fitting  procedure  and  the  traditional  approach,  Figure  13  coverts  Figure  12  to  a  probability- 
square-root  plane.  The  chosen  root  nonnal  parameters  appear  to  provide  an  equally  good  fit  in 
either  case.  Notice  that  &  is  the  cumulative  distribution  function  of  the  standard  normal 
distribution  and  is  its  inverse. 


Figure  13.  Another  view  of  Type  I  vs.  Type  II  root  normal  distributions  compared  to  experimental  results 
obtained  by  Houghton  (1941)  and  Tate  and  Olson  (1962).  The  experimental  results  are  as  reported  by 
Pimentel  et.  al.  (2010). 

9.6  Experimental  Evidence  for  o  =  0.286 

As  one  example,  Spielbauer  et.  al.  (1989)  sprayed  black  liquor  through  splashplate  nozzles  under 
varying  conditions.  They  found  35  different  root  normal  size  distributions  with  a  minimum  value 
of  (7  -  0.230 ,  a  maximum  value  of  a  =  0.330 ,  and  an  arithmetic  average  value  of  a  =  0.286 . 
Over  20%  of  the  tests  obtained  the  average  value  of  a  -  0.286  almost  exactly.  For  additional 
examples  along  these  same  lines,  see  Adams  et.  al.  (1990)  and  Loebker  &  Empie  (2001).  None 
of  the  data  in  Pimentel  et.  al.  (2010)  is  well-fitted  by  a  Type  I  distribution  with  a  =  0.286  ;  thus 
the  plots  seen  in  earlier  cases  are  omitted  here. 
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10.  Conclusions 


As  summarized  in  Table  7,  a  small  number  of  root  normal  size  distributions,  on  the  order  of  a 
half-dozen,  may  provide  reasonably  good  fits  to  a  wide  variety  of  test  data.  No  pretense  is  made 
that  these  are  always  the  best  possible  fits.  In  particular,  in  many  cases,  root  normal  size 
distributions  provide  obviously  poor  fits  to  the  largest  droplets. 


Table  7.  Recommended  Type  IIA  and  Type  IIB  substitutes  for  the  Type  I  distributions  listed  in  Table  5. 


Original 

Recommended 

Type! 

[A  (*) 

Type  I] 

[A  (*) 

Type  IIB  (**) 

C 7 

a 

o 

a 

<7 

a 

1 

0.081 

1 

0.083 

0.9549 

0.086 

0.9963 

2 

0.11 

1 

0.115 

0.9130 

0.125 

0.9922 

3 

0.17 

1 

0.19 

0.7584 

0.24 

0.9708 

4 

0.20 

1 

0.23 

0.6422 

0.34 

0.9389 

5 

0.238 

1 

0.27 

0.5019 

0.47 

0.8650 

6 

0.286 

1 

0.33 

0.2452 

0.72 

0.5738 

(*)- 

Normalized  by  the  MMD  (**)  Normalized  by  the  CMD. 

Unlike  other  size  distributions,  the  parameters  in  root  nonnal  size  distributions  depend  on  which 
average  size,  such  as  the  MMD  or  CMD,  is  used  for  nonnalization.  As  seen  in  Table  7,  the 
difference  is  small  only  when  o  is  small. 

While  not  always  an  exact  match  to  the  traditional  Type  I  distributions,  the  Type  II  and  IIB 
distributions  given  in  Table  7  avoid  severe  singularities  found  in  the  Type  I  distributions,  and 
ensure  true  self-similarity.  Furthermore,  the  preponderance  of  the  experimental  evidence 
reviewed  here  favors  Type  II  over  Type  I  root  normal  size  distributions. 

As  seen  in  Section  9,  most  previous  work  on  root  normal  size  distributions  relies  on  heavily- 
averaged  test  data.  First,  the  results  for  each  test  are  usually  averaged  over  time  and/or  space.  In 
those  rare  cases  where  test  results  are  given  for  specific  points  in  time  and/or  space,  the  results 
often  deviate  substantially  from  the  average,  e.g.,  Ruff  et.  al.  (1992),  Chou  et.  al.  (1997).  Next, 
the  results  are  usually  averaged  across  multiple  test,  sometimes  dozens  or  even  hundreds  of 
different  tests  with  different  outcomes.  When  results  are  given  for  one  specific  test,  the  results 
often  deviate  substantially  from  the  average. 

As  a  typical  conclusion  based  on  such  averages-of-averages,  Spielbauer  et.  al.  (1989)  says  that 
the  “size  distribution  doesn’t  change,  or  changes  very  little,  as  a  function  of  nozzle  geometry, 
flow  conditions,  and  fluid  parameters  . . .  What  is  most  surprising  ...  is  the  similarity  of  the 
droplet  size  distribution  from  the  splashplate  and  swirl  cone  nozzles.”  While  ±15%  or  ±20% 
variations  in  root  nonnal  parameters  such  as  a  may  seem  small,  the  results  given  here  show  that 
they  can  be  highly  significant,  although  this  obviously  depends  on  the  application. 

By  focusing  on  individual  tests  rather  than  ensemble  averages  of  different  tests,  this  treatment 
has  avoided  one  kind  of  averaging.  However,  it  is  hard  to  avoid  the  other  kind  of  averaging  - 
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individual  test  results  may  still  depend  on  averages  in  time  and/or  space.  Accurate  pre-test 
prediction  is  difficult  to  the  extent  that  test  results  depend  on  unknowns  such  as  exactly  when  or 
where  measurements  will  be  taken. 

The  results  presented  here  indicate  that  the  well-known  Simmons  universal  root  normal  size 
distribution  may  not  always  be  especially  common.  In  particular,  for  the  data  compiled  by 
Pimentel  et.  al.  (2010),  Simmons  universal  root  normal  provides  the  best  fit  just  about  10%  of 
the  time;  see  Table  6.  Obviously,  this  observation  only  applies  to  one  particular  collection  of  test 
data.  Different  collections  of  test  data  have,  in  the  past,  yielded  very  different  conclusions. 

Future  work  will  hopefully  provide  additional  insight  into  the  correctness,  completeness,  and 
frequency  of  occurrence  of  the  root  nonnal  parameters  surveyed  here. 

Acknowledgments 

This  work  was  supported  by  the  Nuclear  Effects  Office  of  the  Defense  Threat  Reduction  Agency 
(DTRA).  The  author  wishes  to  thank  Dr.  John  Cockayne  of  Leidos  Inc.  for  his  contributions. 

References 

C.  Aalburg,  B.  Van  Leer,  G.M.  Faeth,  and  K.A.  Sallam,  Properties  of  Nonturbulent  Round 
Liquid  Jets  in  Uniform  Gaseous  Crossflows,”  Atomization  and  Sprays,  15(3),  271-294,  2005 

T.  N.  Adams,  H.  L.  Empie,  N.  Obuskovic,  and  T.  M.  Spielbauer,  Kraft  Black  Liquor  Delivery 
Systems  Report  No.  1,  Institute  of  Paper  Science  and  Technology  Technical  Report,  February 
1990 

N.  Ashgriz,  Editor,  Handbook  of  Atomization  and  Sprays:  Theory  and  Applications,  Springer, 
2011 

E.  Babinsky  and  P.  E.  Sojka,  Modeling  Drop  Size  Distributions,  Progress  in  Energy  and 
Combustion  Science,  28,  303-329,  2002 

L.  P.  Bayvel  and  Z.  Orzechowski,  Liquid  Atomization,  Taylor  and  Francis,  1993 

J.  C.  Bhatia,  J.  Domnick,  F.  Durst,  and  C.  Tropea,  Phase-Doppler  Anemometry  and  the  Log- 
Hyperbolic  Distribution  Applied  to  Liquid  Sprays,  Particle  and  Particle  Systems 
Characterization,  5(4),  153-164,  1988 

C.  E.  Brennan,  Fundamentals  of  Multiphase  Flows,  Cambridge  University  Press,  2005 

W.  K.  Brown,  A  Theory  of  Sequential  Fragmentation  and  Its  Astronomical  Applications,  Journal 
of  Astrophysics  and  Astronomy,  10(1),  89-112,  1989 

W.  K.  Brown  and  K.  H.  Wohletz,  Derivation  of  the  Weibull  Distribution  Based  on  Physical 
Principles  and  Its  Connection  to  the  Rosin-Rammler  and  Lognormal  Distributions,  Journal  of 
Applied  Physics,  78(4),  2758-2763,  1995 


34 


W.-H.  Chou,  L.-P.  Hsiang,  and  G.  M.  Faeth,  Temporal  Properties  of  Drop  Breakup  in  the  Shear 
Breakup  Regime,  International  Journal  of  Multiphase  Flow,  23,  651-669,  1997 

W.-H.  Chou  and  G.  M.  Faeth,  Temporal  Properties  of  Secondary  Drop  Breakup  in  the  Bag 
Breakup  Regime,  International  Journal  of  Multiphase  Flow,  24,  889-912,  1998 

A.  Clauset,  C.  R.  Shalizi,  and  M.  E.  J.  Newman,  Power-Law  Distributions  in  Empirical  Data, 
SIAM  Review,  51(4),  661-703,  2009 

C.  T.  Crowe,  Editor,  Multiphase  Flow  Handbook,  CRC  Press,  2006 

Z.  Dai,  W.-H.  Chou,  and  G.  M.  Faeth,  Drop  Formation  Due  to  Turbulent  Primary  Breakup  at  the 
Free  Surface  of  Plane  Liquid  Wall  Jets,  Physics  of  Fluids,  10,  1 147-1157,  1998 

C.  Dumouchel,  The  Maximum  Entropy  Formalism  and  the  Prediction  of  Liquid  Spray  Drop-Size 
Distribution,  Entropy,  11,  713-747,  2009 

C.  Dumouchel,  J.-B.  Blaisot,  and  V.D.  Ngo,  On  the  Adequacy  Between  the  Laser  Diffraction 
Diameter  Distribution  and  the  3 -parameter  Generalized-Gamma  function,  Chemical  Engineering 
Science,  79,  103,  2012 

H.  J.  Empie,  S.  J.  Lien,  R.  S.  Rumsey,  and  D.  G.  Sachs,  Kraft  Black  Liquor  Delivery  Systems 
Report  No.  4,  Institute  of  Paper  Science  and  Technology  Technical  Report,  January  1993 

H.  J.  Empie,  S.  J.  Lien,  W.  Yang,  and  T.  N.  Adams,  Spraying  Characteristics  of  Commercial 
Black  Liquor  Nozzles,  TAPPI  Journal,  78(1),  121-128,  1995 

H.  J.  Empie,  S.  J.  Lien,  and  W.  Yang,  Drop  Size  Modification  in  Black  Liquor  Sprays  From 
Commercial  Nozzles  Using  Pulsed  Flow,  Atomization  and  Sprays,  7,  457-466,  1997 

S.  K.  Friedlander  and  C.  Wang,  The  Self-Preserving  Particle  Size  Distribution  for 
Coagulation  by  Brownian  Motion,  Journal  of  Colloid  and  Interface  Science,  22,  126-32,  1966 

H.  G.  Houghton  in  J.  H.  Perry,  Editor,  Perry ’s  Chemical  Engineers  ’  Handbook,  2nd  ed. , 
McGraw-Hill,  1941 

L.-P.  Hsiang  and  G.  M.  Faeth,  Near-Limit  Drop  Deformation  and  Secondary  Breakup, 
International  Journal  of  Multiphase  Flow,  18,  635-652,  1992 

L.-P.  Hsiang  and  G.  M.  Faeth,  Drop  Deformation  and  Breakup  Due  to  Shock  Wave  and  Steady 
Disturbances,  International  Journal  of  Multiphase  Flow,  21,  545-560,  1993 

D.  D.  Joseph,  J.  Belanger,  and  G.  S.  Beavers,  Breakup  of  a  Liquid  Drop  Suddenly  Exposed  to  a 
High-Speed  Airstream,  International  Journal  of  Multiphase  Flow,  25,  1263-1303,  1999 


35 


D.  D.  Joseph,  G.  S.  Beavers,  and  T.  Funada,  Rayleigh-Taylor  Instability  of  Viscoelastic  Drops  at 
High  Weber  Numbers,  Journal  of  Fluid  Mechanics,  453,  109-132,  2002 

K.  Lee,  C.  Aalburg,  F.  J.  Diez,  G.  M.  Faeth,  and  K.  A.  Sallam,  Primary  Breakup  of  Turbulent 
Round  Liquid  Jets  in  Unifonn  Crossflows,  AIAA  Journal,  45(8),  2007 

K.  W.  Lee,  H.  Chen,  and  J.  A.  Gieske,  Log-Normally  Preserving  Size  Distribution  for  Brownian 
Coagulation  in  the  Free-Molecular  Regime,  Aerosol  Science  and  Technology,  3,  53-62,  1984 

K.  E.  J.  Lehtinen  and  M.  R.  Zachariah,  Self-Preserving  Theory  for  the  Volume  Distribution  of 
Particles  Undergoing  Brownian  Coagulation,  Journal  of  Colloid  and  Interface  Science,  242,  3 14- 
318,2001 

Li,  X.  and  Tankin,  R.  S.,  Droplet  Size  Distribution:  A  Derivation  of  a  Nukiyama-Tanasawa  Type 
Distribution  Function,  Combustion  Science  and  Technology,  65-76,  1987 

H.  Liu,  Science  and  Engineering  of  Droplets:  Fundamentals  and  Applications,  William  Andrew 
Publishing,  2000 

D.  W.  Loebker  and  H.  J.  Empie,  Independently  Controlled  Drop  Size  in  Black  Liquor  Sprays 
into  the  Kraft  Recovery  Boiler  Using  Effervescent  Atomization,  Journal  of  Pulp  and  Paper- 
Science,  27(1),  18-25,  2001 

B.  Miller,  K.A.  Sallam,  M.  Bingabr,  K.-C.  Lin,  and  C.  Carter,  Breakup  of  Aerated  Liquid  Jets  in 
Subsonic  Crossflow,"  Journal  of  Propulsion  and  Power,  24(2),  253-258,  2008 

R.  Pimentel,  R.  Stowe,  P.  Harris,  A.  DeChamplain,  and  D.  Kretschmer,  Spray  Characterization 
Based  on  the  Pearson  System  of  Frequency  Curves,  Atomization  and  Spravs,  20(5),  365-386, 
2010 

T.  Paloposki  and  N.  E.  Fagerholm,  The  Atomization  of  Heavy  Fuel  Oil  with  Swirl  Pressure  Jet 
Atomizer,  Helsinki  University  of  Technology,  Institute  of  Energy  Engineering  Research  Report 
TKK-KO/ET-13,  NTIS  Accession  Number  DE87751 107,  1986 

T.  Paloposki,  Drop  Size  Distributions  in  Liquid  Sprays,  Acta  Polytechnica  Scandinavica, 
Mechanical  Engineering  Series  114,  Helsinki,  Finland,  1994 

G.  A.  Ruff,  P.-K.  Wu,  L.  P.  Bernal  and  G.  M.  Faeth,  Continuous-  and  Dispersed-Phase  Structure 
of  Dense  Nonevaporating  Pressure- Atomized  Sprays,  Journal  of  Propulsion  and  Sprays,  8(2), 
1992 

K.  A.  Sallam,  Z.  Dai,  and  G.  M.  Faeth,  Drop  Formation  at  the  Surface  of  Plane  Turbulent  Liquid 
Jets  in  Still  Gases,  International  Journal  of  Multiphase  Flow,  25,  1161-1180,  1999 


36 


K.A.  Sallam,  C.  Aalburg,  G.M.  Faeth,  K.-C.  Lin,  C.D.  Carter,  and  T.A.  Jackson,  Primary 
Breakup  of  Round  Aerated-Liquid  Jets  in  Supersonic  Crossflows,  Atomization  and  Sprays,  16(6), 
657-672,  2006 

H.  C.  Simmons,  The  Correlation  of  Drop-Size  Distributions  in  Fuel  Nozzle  Spray,  Parts  I  and  II, 
Journal  of  Engineering  for  Power,  99,  309-319,  1977 

P.  T.  Spicer  and  S.  E.  Pratsinis,  Coagulation  and  Fragmentation:  Universal  Steady-State  Particle- 
Size  Distribution,  AIChE  Journal,  42,  1612-1620,  1996 

T.  M.  Spielbauer,  T.  N.  Adams,  J.  E.  Monacelli,  and  R.  T.  Bailey,  Droplet  Size  Distribution  of 
Black  Liquor  Sprays,  1989  International  Chemical  Recovery  Conference,  Ottawa,  Canada, 
TAPPI/CPPA,  1989 

R.  W.  Tate  and  W.  R.  Marshall  Jr.,  Atomization  by  Centrifugal  Pressure  Nozzles,  Chemical 
Engineering  Progress,  49,  169-174,  1953 

R.  W.  Tate  and  E.  O.  Olson,  Spray  Droplet  Size  of  Pressure- Atomizing  Burner  Nozzles, 

ASHRAE  Journal,  4(3),  39-43,  1962 

J.  M.  Tishkoff,  GMR  Laser-video  Imaging  System  Results  for  the  ASTM  Round- 
Robin’  Nozzle  Test,  General  Motors  Research  Publication  GMR-3098,  1979 

G.  M.  Turner  and  R.  W.  Moulton,  Drop  Size  Distributions  from  Spray  Nozzles, 

Chemical  Engineering  Progress,  49(4),  185-190,  1953 

T.  G.  Theofanous,  Aerobreakup  of  Newtonian  and  Viscoelastic  Liquids,  Annual  Review  Fluid 
Mechanics,  43,  661-690,  2011 

P.-K.  Wu,  G.  A.  Ruff,  and  G.  M.  Faeth,  Primary  Breakup  in  Liquid-Gas  Mixing  Layers, 
Atomization  and  Sprays,  1,  421-440,  1991 

P.-K.  Wu,  L.-K.  Tseng,  and  G.  M.  Faeth,  Primary  Breakup  in  Gas/Liquid  Mixing  Layers  for 
Turbulent  Liquids,  Atomization  and  Sprays,  295-317,  1992 

P.-K.  Wu  and  G.  M.  Faeth,  Aerodynamic  Effects  on  Primary  Breakup  of  Turbulent  Liquids, 
Atomization  and  Sprays,  3,  265-289,  1993 

P.-K.  Wu,  R.  F.  Miranda,  and  G.  M.  Faeth,  Effects  of  Initial  Flow  Conditions  on  Primary 
Breakup  of  Nonturbulent  and  Turbulent  Round  Liquid  Jets,  Atomization  and  Sprays,  5,  175-196, 
1995 


37 


DISTRIBUTION  LIST 
DTRA-TR-16-3 


DEPARTMENT  OF  DEFENSE 


DEFENSE  THREAT  REDUCTION 
AGENCY 

8725  JOHN  J.  KINGMAN  ROAD 
STOP  6201 

FORT  BELVOIR,  VA  22060 
ATTN:  CAPT  PAUL 
CLEMENT 

DEFENSE  TECHNICAL 
INFORMATION  CENTER 
8725  JOHN  J.  KINGMAN  ROAD, 
SUITE  0944 

FT.  BELVOIR,  VA  22060-6201 
ATTN:  DTIC/OCA 

DEPARTMENT  OF  DEFENSE 
CONTRACTORS 

QUANTERION  SOLUTIONS,  INC. 
1680  TEXAS  STREET,  SE 
KIRTLAND  AFB,  NM  87117-5669 
ATTN:  DTRIAC 


DL-1 


